1 



On 
OS 
t— i 



Calculation of Moments of Nucleon Structure Functions* 

R.C. Brower, S. Huang, J.W. Negele, A. Pochinsky and B. Schreiber a 

a Center for Theoretical Physics, LNS, Massachusetts Institute of Technology, 
Cambridge MA 02139, USA 

Preliminary results are presented in our program to calculate low moments of structure functions for the proton 
and neutron on a 24 3 x 32 lattice at /3 = 6.2. A comparison is made for a variety of smeared nucleon sources and 
preliminary results for the calculation of the nucleon tensor charge are presented. 



> 

On 

\o 

O 

OO 
O 

On 



X 



1. INTRODUCTION 

There are now extensive experimental results 
on nucleon structure functions in deep inelastic 
scattering, providing a detailed empirical knowl- 
edge of the distribution of quarks and gluons 
in the nucleon. Thus a major theoretical chal- 
lenge for lattice QCD is to understand this data 
from first principles and assess errors in current 
methods due to the quenched approximation and 
the finite lattice spacing in naive lattice opera- 
tors. We envision a gradual improvement of lat- 
tice methodology to reduce these errors. Since 
structure functions cannot be calculated directly 
on a Euclidean lattice, we consider low moments 
which are related to matrix elements of local op- 
erators by the operator product expansion. The 
lattice matrix elements relevant to moments of 
Pi, gi and hi structure functions are being cal- 
culated on quenched 24 3 x 32 lattices at p = 6.2 
at three quark masses, and initial results are pre- 
sented below. 



2. GENERAL FORMALISM ON THE 
LATTICE 

In the continuum, the quark contribution to 
the moments of the unpolarized Pi, longitudi- 
nally polarized g\ and transversely polarized h\ 
structure functions are related to the matrix ele- 
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ments, 

S 

(pstywsD^ . . . Df, n tp\ps) ~ a n (n), 
(ps^dv^D^ . . . D^ijj\ps) ~ t n (p), 

respectively, where we have indicated on the right 
the scalar coefficient in the notation of Refs. M. 

Our methodology for calculating these opera- 
tors on the lattice is similar to Gockeler et aljjj. 
The covariant derivative is replaced by a lattice 
finite difference and the operators are classified by 
representations of the surviving discrete subgroup 
Hi of the continuum Euclidean Lorentz group 
50(4) The renormalization factors Zq are 
calculated using a new method j^] to one loop or- 
der in lattice perturbation series to account for a 
multiplicative shift of the Wilson coefficient due 
replacing the MS scheme at scale /x by the lattice 
spacing a. Our independently derived results con- 
firm those in the literature and extend them 
where necessary. For our current set of moments 
the results are given in the table. 

Typical elements of an irreducible represen- 
tation (IR) of Hi are given in the second col- 
umn, where the {■ • •} and [■ ■ •] brackets imply 
symmetrized and anti-symmetries indices respec- 
tively. In many cases there is more than one IR 
corresponding to the continuum operator. We are 
exploring other representations to have a redun- 
dant set of measurements. The renormalization 
constants defined by the formula, 



Za = 1 
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C F [ 7G ln(^a) + Bq] 
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are given in the third and the fourth column. 
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3. PRELIMINARY RESULTS 

On the CM-5 at MIT, we have generated 150 
independent 24 3 x 32 lattices using the stan- 
dard Wilson action. Dirac propagators for Wil- 
son fermions with r = 1 are being calcu- 
lated by conjugate gradient iterations with red- 
black preconditioning for the values of K = 
0.15200,0.15246,0.15294 corresponding to m, « 
150, 98 and 45 MeV. Our estimates for k c and the 
lattice spacing derived from 7r and p masses are 
in a good agreement with those of UKQCD. 

3.1. Hadronic Sources. 

To compute the matrix elements we consid- 
ered four types of sources: point sources (P), 
gauge invariant Wuppertal || sources (W) and 
both Coulomb gauge fixed Wuppertal (U) and 
Gaussian (G) smeared sources. To determine 
the most suitable form of the source, we inves- 
tigated the plateau in \a.{G{t)/G(t + 1)) for the 
two point functions for the pion, rho and nu- 
cleoli sources. For the proton creation operator, 
we used J^x) = e abc u^(x)u b a (x)T a ^d c p(x), with 
r = C75 and with all quark operators truncated 
to 2 upper components. 

Some comparisons are shown for the nucleon 
case in Fig. |l|. The Gaussian (G) and the two 
Wuppertal (W & U) sources were adjusted so 



yj (x 2 ) ~ 6.7a for each quark field, since this 
smearing produced the least noisy results in all 
three cases. For more localized sources the ex- 
cited states are more prominent, whereas for less 
localized sources the signal becomes noisier at 
large distances. 

As seen in Fig. |], smearing both the sources 
and the sink results in substantially noisier behav- 
ior than smearing only the source. On the scale of 
the errors in Fig. [ll there is no significant differ- 
ence between smeared source-point sink vs. point 
source-smeared sink. It is interesting to note that 
the gauge fixed Wuppertal (D) and gauge invari- 
ant Wuppertal (W) sources are essentially equiv- 
alent. 
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Figure 1. Effective mass as a function of distance 
for the point-point (PP), Wuppertal-point (WP), 
gauge fixed Wuppertal-point (UP), Gaussian- 
Gaussian (GG) and Gaussian-point (GP) nucleon 
correlators. 

3.2. The Tensor Charge 

We illustrate our preliminary results with 
a first (low statistics) measurement of the 
tensor charge M. (Also see results by 
Aoki et al@.) To date we have analyzed 
15,10,5 configurations with p = (0,0,0) for 
k = 0.15200,0.15246,0.15294 respectively and 
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20,15,15 configurations at p = (1,1,0) for the 
same kappa values. The Gaussian smeared source 
is placed at t = 8 and the momentum projected 
point sink at t\ = 24. 

We measure three point functions in the fol- 
lowing way. Once a gauge field is generated, it 
is gauge fixed to the Coulomb gauge and a set 
of quark propagators is computed. Then for ev- 
ery operator the corresponding 3-point function 
is constructed from the set of propagators. We 
calculate separately the 2-point function, compar- 
ing backward and forward propagators indepen- 
dently to verify convergence. Finally, we use the 
jackknifc method to estimate statistical errors of 
the ratio of three- and two-point functions. Since 
we need to fix the final state momentum of the 
nucleon before calculating the backward propa- 
gators, the set of propagators has a built-in final 
state momentum. 

The first moment t\ of the transversity distri- 
bution h\ (x) can be measured using the final state 
at rest p = (0,0,0). In Fig. || the three upper 
plots show preliminary results for t\ at the three 
kappas we used. Even for the low statistics (15 
configurations for the smallest k), the plateau in 
the signal is clearly visible. One can also see that 
as k approaches the chiral limit, the signal be- 
comes noisier as expected. 

Since other moments require non-zero momen- 
tum at the sink, it might be useful to consider a 
single set of propagators with non-zero momenta 
to measure all observables. As a comparison, 
the lowest plot in Fig. ||D gives the same ten- 
sor charge computed at p = (1,1,0). To make 
this a fair comparison with the p = (0, 0, 0) case, 
we have used the same set of configurations with 
propagators calculated to the same precision. It 
appears to have considerably stronger fluctua- 
tions, thus raising the question as to whether the 
exclusive use of non-zero momenta is an optimal 
strategy. 
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Figure 2. u and d contributions to t\ for p = 
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